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Abstract 



The recent reported charmed-strange resonance at 2.32 GeV/c suggests 
a possible multiquark state. Three types of multiquark bound states are 

m 
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reviewed. A previous model-independent variational approach considers a 
tetraquark with two heavy antiquarks and two light quarks as a heavy antidi- 
quark with the color field of a quark bound to the two light quarks with a wave 
function like that of a heavy baryon. Results indicate that a charmed-strange 
D ! tetraquark csud or a bottom-strange tetraquark bsud with this "diquark- 

heavy-baryion" wave function is not bound, in contrast to "molecular-type" 
D — K and B — K wave functions. However, a charmed-bottom tetraquark 
^ ■ cbud might be bound with a very narrow weak decay mode. A "molecular- 

type" D — B state can have an interesting B c ir decay with a high energy 
pion, 



I. THREE TYPES OF MULTIQUARK STATES 

The recent observation [1] of a charmed-strange state at 2.32 GeV that decays into D s ii 
suggests a possible four-quark state (tetraquark). [2-8] 
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Three different mass scales are relevant to the description of multiquark hadrons, the 
nuclear-molecular scale, the hyperfine or color-magnetic scale and the diquark scale. 

The nuclear scale is characterized by the deuteron, a bound state of two color singlet 
hadrons with a reduced mass of 500 MeV, a binding energy of several MeV and a radius of 
ps M w . The underlying quark structure of the hadrons plays no role. The kinetic energy of 
the state confined to this radius is 

T N = p 2 /M N ps Ml/M N Ps 20MeV (1.1) 

No two-meson bound state containing a pion has been found. The reduced mass of any 
such state is too small to be bound in a radius of ps M n by a similar interaction; its kinetic 
energy would be too high. 

T w Ps Ml /M w ps 140 MeV (1.2) 
The two-kaon system with a reduced mass of 250 MeV seems to be on the borderline, 

T K = p 2 /M K ps Ml/M K ps 40 MeV (1.3) 

Suggestions that the f and a Q mesons are deuteron-like KK states or molecules are 
interesting, but controversial. There is no unambiguous signature because KK couples to 
n — n and i] — n and both states break up strongly. 

The D — K system with a kinetic energy 

T DK = p 2 (M D + M K )/2M D M K ps Ml{M D + M K )/2M D M K ps 25 MeV (1.4) 

is therefore an attractive candidate for such a state [3-8] . The transition for the 1 = DK 
state to Dstt is isospin forbidden; thereby suggesting a narrow width. 

The color-magnetic scale is characterized by a mass splitting of the order of 400 MeV; 
e.g. the K* — K splitting. Recoupling the colors and spins of a system of two color-singlet 
hadrons has been shown to produce a gain in color-magnetic energy [2-4]. However, whether 
this gain in potential energy is sufficient to overcome the added kinetic energy required for 
a bound state is not clear without a specific model. 
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The diquark scale arises when two quarks are sufficiently heavy to be bound in the well 
of the coulomb-like short-range potential required by QCD. A heavy antidiquark in a triplet 
of color SU(3) has the color field of a quark and can be bound to two light quarks with 
a wave function like that of a heavy baryon. Since the binding energy of two particles in 
a coulomb field is proportional to their reduced mass and all other interactions are mass 
independent, this diquark binding must become dominant at sufficiently high quark masses. 

II. THE DIQUARK-HEAVY-BARYON MODEL FOR TETRAQUARKS 

We now examine the diquark-heavy-baryon model for states containing heavy quarks. 
Our "model-independent" approach assumes that nature has already solved the problem of 
a heavy color triplet interacting with two light quarks and given us the answers; namely 
the experimental masses of the A, A c and A 6 . These answers provided by nature can now 
be used without understanding the details of the underlying theoretical QCD model. This 
approach was first used by Sakharov and Zeldovich [9] and has been successfully extended 
to heavy flavors [10-12]. 

The calculated mass can be interpreted as obtained from a variational principle with 
a particular form of trial wave function [5]. This model neglects the color- magnetic inter- 
actions of the heavy quarks, important for the charmed-strange four-quark system at the 
colormagnetic scale [2-4] and is expected to overestimate the mass of a csud state. Thus 
obtaining a model mass value above the relevant threshold shows only that this type of 
diquark-heavy-baryon wave function does not produce a bound state; i.e that the heavy 
quark masses are not at the diquark scale. The previous results [3,4] at the colormagnetic or 
nuclear-molecular scale should be better. However, the be system may already be sufficiently 
massive to lead to stable diquarks and the model predictions for the cbud state may suggest 
binding. 

We first apply this model to a csud state with a light ud pair seeing the color field of 
the cs antidiquark like the field of a heavy quark in a heavy baryon. The cs antidiquark 
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differs from the cs in the D s by having a QQ potential which QCD color algebra requires [5] 
to have half the strength of the QQ potential in the D s . The tetraquark mass is estimated 
by using the known experimental masses of the heavy baryons and heavy meson with the 
same flavors and introducing corrections for the difference between the heavy meson and the 
heavy diquark. 

M(csud) = m c + m s + m u + m d + (H udQ ) + (H ud ) + (T cs ) cs + (V cs ) cs (2.1) 
M(cs) =m c + m s + (T cs ) cs + (V cs ) cs (2.2) 
M(D S ) = m c + m s + (T cs ) cS + (V cS ) cS (2.3) 
M (A) =m s + m u + m d + (H ud ) + (H udQ ) (2.4) 

M(A C ) =m c + m u + m d + (H ud ) + (H udQ ) (2.5) 

where H ud and H ud g respectively denote the Hamiltonians describing the internal motions 
of the ud pair and of the three-body system of the ud pair and the antidiquark which 
behaves like a heavy quark, T cs and V cs denote the kinetic and potential energy operators 
for the internal motion of a cs diquark which is the same as that for a cs antidiquark. 
The expectation values are taken with the "exact" wave function for the model, with the 
subscript cs indicating that it is taken with the wave function of a diquark and not of the 
D s . The kinetic energy operator T cs is the same for the cs diquark and the D s but the 
potential energy operators V cs and V cS = 2V CS differ by the QCD factor 2. This difference 
between cs diquark and D s wave functions is crucial to our analysis. 

The quark masses m q are effective constituent quark masses and not current quark 
masses. We follow the approach begun by Sakharov and Zeldovich [9] who noted that both 
the difference m s — m u between the effective masses of strange and nonstrange quarks and 
their ratio m s /m u have the same values when calculated from baryon masses and meson 
masses. 
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(m s - m u ) Bar = M A -M N = 177 MeV (2.6) 



/ \ 3(M K * - M p ) + M K - M n 18nMV (07 , 
{m s - m u ) Me s = —. = 180 MeV (2.7) 



M -- M « L5 3 = M '~ M ' =1.61 (2.8) 



^uJ Bar M s * -M s \m u J Mes M K * - M K 

where the "Bar" and "Mes" subscripts denote values obtained from baryons and mesons, 
respectively. Similar results have since been found for hadrons containing heavy quarks 
along with many more relations using these same effective quark mass values for baryon 
magnetic moments and hadron hyperfine splittings [10-12]. We therefore assume that the 
values of the effective quark masses m q remain the same for all meson and baryon states in 
our analysis. 

Substituting eqs.(2.2 - 2.5) into eq. (2.1) gives 

M{csud) = (1/2) • [M(D S ) + M(A) + M(A C )] + (5H CS ) (2.9) 

where (SH CS ) expresses the difference between the D s and the cs wave functions 

(5H CS ) = (T cs ) cs + (V cs ) cs - (1/2) • [(T cs ) c - S + (V cS ) cS ] (2.10) 

To calculate (5H CS ) we first improve on the treatment of ref [5] and define the Hamiltonian 

H(a) = aT cs + V cs = aT cs + (1/2) • V cS (2.11) 

This Hamiltonian H(a) is seen to describe both the cs diquark and the D s 

M(cs)=m c + m s + (H(a)) a=1 (2.12) 

M(D S ) =m c + m s + 2- (if (a)) Q =(i/ 2 ) (2.13) 



(5H CS ) = (H(a)) a=1 - (H(a)) a={1/2) (2.14) 

To evaluate (SH CS ) we use the Feynman-Hellmann theorem and the virial theorem to obtain, 
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q: 



This expression can be simplified by using the Quigg-Rosner logarithmic potential [13] with 
its parameter V Q determined by fitting the charmonium spectrum. 

V c Q s R = (1/2) • V ■ log(r/r ) (2.17) 

(5H cs ) QR =^f — = ^2 = 126 MeV (2.18) 
Substituting experimental values then shows M(csud) well above the -D.K" threshold [14] 
M(cswd) = (1/2) • [M(D S ) + M (A) + M(A C )] + (6# cs ) = 2685 + 126 = 2811 MeV (2.19) 

M(csud) = 2811 MeV > M(D) + M{K) = 2361 MeV (2.20) 

In the limit of very high heavy quark masses this model must give a stable bound state. 
The cs diquark is evidently not heavy enough to produce a bound diquark-heavy-baryon 
state. 

A similar calculation for bsud indicates that the 6s diquark is also not heavy enough. 
M(bsud) = (1/2) • [M(B S ) + M (A) + M(A b )} + (5H bs ) = 6180 MeV > M(B) + M(K) = 5773 MeV 

(2.21) 

However, the be diquark may be heavy enough to produce a bound four-quark state. 
M(cbud) = (1/2) • [M(B C ) + M(A b ) + M(A C )] + (5H CS ) = 7280 ± 200 MeV (2.22) 

M(D) + M(B) = 7146 MeV (2.23) 

Here the experimental error on the B c mass is too large to enable any conclusions to be 
drawn. But if the bound state exists, it may produce striking experimental signatures. 
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A bound cbud, cbuu or cbdd state would decay only weakly, either by b-quark decay 
into two charmed mesons (with the same sign of charm, so that there cannot be a J/psi 
decay mode), or a c-quark decay into a b meson and a strange meson. The signature with 
a vertex detector will see a secondary vertex with a multiparticle decay and one or two 
subsequent heavy quark decays and either one track or no track from the primary vertex to 
the secondary. 

On the other hand, if the 2.32 GeV state seen by BaBar is really a DK 1 = molecule 
with an isospin violating D s — n decay, the analog for the be system is a BD molecule with 
either / = 1 or I — and a B c — it decay, which is isospin conserving for / = 1 or isospin 
violating for I = 0. 

Here the masses are very different and give a completely different signature with a high 
energy pion. M(B) = 5279 MeV, M(D) = 1867 MeV. This gives M(B) + M(D) = 7146 
MeV, while M(B C ) = 6400 ± 400 MeV. So a molecule just below BD threshold would just 
rearrange the four quarks into Bc-pi and fall apart, either with or without isospin violation, 
giving a neutral or charged pion having a well defined energy of 750 ± 400 MeV with the 
precision improved by better measurements. 

In any case this is a striking signal which cannot be confused with a qq state. Experiments 
can look for a resonance with a pion accompanying any of B c states. 

It is a pleasure to thank E. L. Berger, T. Barnes, F. E. Close, M. Karliner, J. Napolitano 
and V. Papadimitriou for helpful discussions. 
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